NOTES FROM THE SEARCH
9 kids

Population with 50 pre-zygotic and 100 post zygotic as above seems stable at 200.

On run 4 same but only 8 kids is stable.  Saved for now as InPop. 

6 kids:  200/4000 pops, who lived?

Running at ID, then 1 sim, 1,000 gens, 6 kids, 0 init pop, varying max, 0 genes, 0 gene mutations, 50 pre-zygotic sites, 2,000 epimuts, 23 impact, 100 post-zygotic sites, s,000 epimuts, 23 impact, 1 delay, add file “unknown” or “unkn,” varying number.  Max pop across, imported pop down.  How many lived out of 10? Accumulating 6, work sheet 1.

Second go at 1300 and 35 with 6 offspring showed no inbreeding depression, sheet 1 run 18.  Need to reduce preZ more.  Also will increase kids from 6 to 10.  Impact stays at 35.  Tweak epimut rate starting at 2500.  Drop pre-z sites from 50 to 25.  ID 19, sheet 2, 10 sims, 1k gens, pop init and top 200, 0 genes and mutations, 1 delay as always thus:   lived at 200/lived at 4k. 

10 sims, 1k gens, 11 kids, 0 genes and mut rate, 12 pre-z sites, listed pre-z change rate and impact, 100 post-z sites, same change and impact, 1 delay as ever, “!” = added file.

Now we go one with the same numbers but vary max pop to 1k, 2k, 3k.  1k worked twice J 
Let’s wait on that.

10 sims, 1k gens, 12 kids, 100 init popsize, 100 max popsize, 0 genes, 0 muts, 20 preZ sites, 5k change, 50 impact, 100 postZsites, 5k change, 50 impact, 1 delay, live at 100  / live at 1000, prepped 100 as run 45.

2 peaks looked just a little better (slightly slower climb and slightly more complex fall) than run 18 acc 6, which had pre-z sites 12 and kids 11

I shall double the inpop size by adding the same numbers twice.  ;o). 

Hmmm. Not so good. 

All right, the grand mix of acc 7 run 6 was a disaster.  So I’ve pulled up the parameters for “double Run” back in 2012.  They go:
200 kids, 100 pre sites, 12,000 change rate, 1,200 impact and 100 post sites, 1,200 change rate and 1,200 impact.  If I divide them all by 4, that ought to work the same (it really doesn’t work, but I suspect it’s close.) as 50 kids. 

That gives us 24 kids, 3,000 change, 300 impact pre and 300 change and 300 impact post.  We’ll run that for 200/2000.  Acc 7 work sheet 1 run 7. meth rate and impact will be pre/post

So we’ll do a prep run at 200 and then bring them into the 2000 regime. 
Tomorrow. … That would be today, Sep 18, 2016
Did the prep run; ok.  Then tried the main run; dead and. I suspect a problem with the program.  So here we go; I’ll do the prep run (Acc 7, run 12 starting with source code.  (Yeah, sure, I’ve rebooted a lot.)

Here it is; lasted a long time. 
1379 9267 13167 13286 13416 13180 12808 12844 12825 12620 12635 12535 12396 12203 11917 11905 11967 11971 11721 11613 11454 11346 11131 11087 11066 10707 10367 10215 10014 9899 9917 9604 9419 9163 9191 9195 9336 9041 8982 9026 8849 8704 8676 8361 8327 8408 8330 8166 7915 7921 7956 7943 7759 7758 7636 7618 7564 7157 6909 6889 6354 6393 6245 6363 6566 6259 6292 6373 5942 5535 6018 6031 5730 5943 5572 6003 5754 5068 5101 5361 5706 5409 5721 5637 5904 5825 6595 6483 6561 6581 5996 6416 6868 7397 7535 7656 8022 7696 8323 8017 7856 7127 6890 8305 8051 7026 7525 7796 8019 7988 8323 7305 4904 7892 7220 6673 7140 5575 8404 8212 8028 8348 8271 8106 7451 7908 8110 8031 7899 7731 7638 7444 7245 7184 6766 6542 5778 5723 5838 6623 6488 5325 6194 6315 6311 6307 6534 6609 6479 6709 6804 6400 6452 6808 7005 6996 6954 6228 7379 7410 7469 7608 7303 6977 6562 6533 6127 6790 6948 6110 6835 6800 6970 6956 6546 5553 6244 6490 6974 7045 7274 7091 7302 7493 7289 7374 7392 7219 7423 6856 7001 6760 6834 7033 6918 6627 6655 6693 6821 6833 6852 6442 6490 6410 6841 6966 7330 7526 7619 7525 7556 7638 7678 7998 8099 8292 8504 8714 8868 8767 8731 8708 8943 8918 9337 8954 8849 8892 8589 8662 8423 8286 8394 8619 8553 8465 8314 8141 8200 8126 7884 7883 7748 7481 7434 7301 7057 7011 6906 6997 6960 6809 6717 6597 6420 6486 6143 6359 6160 6112 5756 5425 5683 5252 5372 5069 5527 5449 5477 5664 5301 5492 5725 5761 6126 5976 6131 6313 6405 6587 6141 6477 6848 6872 6204 7392 7480 7451 7466 7456 7044 7308 6915 7022 7109 7296 7183 6801 7000 6939 6829 6715 6826 6293 6029 5782 4178 6143 5962 4220 6218 4650 7119 6951 6557 7097 7028 6931 6237 6548 6647 6053 6676 6429 6575 6426 6128 6011 5916 5874 5856 5731 5725 5791 5609 5253 5204 4968 5302 5137 4908 4564 4939 4809 4750 4822 5185 5261 5168 5441 5165 5721 6621 6676 7136 7729 7958 7234 7524 7316 5906 5441 6802 5929 6897 6136 6755 7709 8202 1993 0

Ixnay.  It’s just pre-zygotic.  I’ll makethe Pre to 100 sites and see.  200/2000

500 looks OK so I’ll graph the last simulation as “Last sim Acc 7 run” 20

The graph says we need more pre-Z so we go to pre-Z 600, 600 and leave post-Z 500, 500. … 6 died, so drop the post z to 400, 400. So 5 died; we drop her to 350, 350 and 8 died.  Maybe we should go back to* 400, 400.  

But we’ll go down to 300, 300. 9 died in the first generation.  Not plausible.  Pewter is down.  I’ll shut her down for the night. 

Trying again on Dec. 18, 2016
Last run was (worksheet2.doc run 5): ID 23; Simulations 10; maximum generations 1k; max offspring 25; initial population 200; Max population 200; genes per chromo 0; gene mutations 0; pre-zygotic sites 100, pre-zygotic mutations 600; pre-zygotic impact 600; post-zygotic sites 100; post-zygotic mutations 350; post-zygotic impact 350; delay 1; add file no.  ONLY 2 SURVIVED. But at 500/500, 500/500 starting new population of 200, all lived. 
So we shall try to save a population of 200 with those parameters and max pop 200.  It will be on worksheet 2 column 7.

OK, now it’s Jan 8, 2017, but we’re continuing as planned: ten of ten runs survived and I saved 200. 

We shall run it again, bringing in the saved 200 and raising max pop to 1,000, 10 reps as in column 8 work sheet 2.   OK then 2 died.  We hold everything the same but go up to4,000.

And that will be worksheet 2, column 9.So it was, and they all died, a couple of 10 early at 88 and 89.  So now we just grind away single runs till we get on dying then, maybe about five runs.

Worksheet 2,column 10 will be just the same but 1 rep and we’ll bring in 200 from inpop. Sooo it died out after 5 peaks and about 200 generations; column 11, run 29 on worksheet 2 will be just the same. Alas, all the work with a balky keyboard to make that graph and the computer et it.  L

We’ll do a few more runs like column 10

Run 30 had three rough peaks and lasted 109 generations.

Run 31 did not have clear peaks and lasted 116 generations.

Run 32 did not have good peaks and lasted 180 generations. 

Run 33 lasted 111 generations.  We’ll have a look.  Not too shabby. J

Run 34 lasted 200 generations.  No need to look at it. 

Run 35 lasted 200 generations.  It’s all a hoax, this reality thing, you know?  We are running Accum 5 inside the compiler. 

Run 36 survived. Wtf
Run 37 died at generation 171.

Run 38, 39 and 40 all lived.  L
Neither 41 nor 42 failed early, so I ran 43 with 200 new population and 200 seasoned population.  It collapsed right away, as did 10 runs in 44, although it looks like it might be a bug.  45 will be 200 fresh ones and 100 seasoned ones.  … exactly the same: first sim died after 2 generations and the rest died after one.  Sooo … we drop the 100 t0 64 in 46 and grope from there.  Nope; on to 32.

Run 47 column 10 worksheet 3 5 died, last at gen 73. Let’s look. OK, seems like it worked. J

Plunging ever onward we’ll pull in 200 of an earlier run that had survived for 1,000 generations and run it for another 200 gens as run 48 on worksheet 3 column 11.  It lived.  Let’s look.  It oscillates in a steady way.  Good enough. J

All right.  Three good and one to go.  Need damped oscillation.  Run 49, which will be worksheet 3, column 12, will be just like 48 but start with 100 pulled in from the old population.  Good.  Now the damped oscillation begins to appear.  J

(what comes through as "J" is pretty much a smiley face)

Home page.